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Abstract. We assess the theoretical uncertainties on the total charm cross section. 
We discuss the importance of the quark mass, the scale choice and the parton 
densities on the estimate of the uncertainty. We conclude that due to the small 
charm quark mass, which amplifies the effect of varying the other parameters in 
the calculation, the uncertainty on the total charm cross section is difficult to 
quantify. 

1 Introduction 

Open charm measurements date back to the late 1970s when D and D mesons were first de- 
tected, completing the picture of the fourth quark begun when the J/ip was detected in pBc 
and e + e~ interactions. The charm quark was postulated to have a mass between 1.2 and 1.8 
GeV, within the regime of perturbative quantum chromodynamics (pQCD). Because of its 
rather large mass relative to the u, d and s quarks, it is possible to calculate a total cc cross 
section, not the case for lighter flavors. Charm hadrons are usually detected two ways. The 
reconstruction of decays to charged hadrons such as D° — > K~tt + (3.8%) and D + — > K^tt + tt + 
(9.1%) gives the full momentum of the initial D meson, yielding the best direct measurement. 
Charm can also be detected indirectly via semi-leptonic decays to leptons such as D — > Klvi 
although the momentum of the parent D meson remains unknown. Early measurements of 
prompt leptons in beam dump experiments assumed that the density of the dump was high 
enough to absorb semi-leptonic decays of non-charm hadrons, leaving only the charm compo- 
nent. At modern colliders, it is not possible to use beam dumps to measure charm from leptons 
but, at sufficiently high pt, electrons from charm emerge from hadronic cocktails |H2j . 

Although D mesons alone are often used to calculate the total cc cross section, other charm 
hadrons also exist. The excited D states, D*s, decay primarily to charged and neutral D mesons. 
The charm-strange meson, the D s , decays to charged hadrons as well as semi-leptonically. The 
lowest mass charm baryon is the A£ which decays primarily to A(uds) but also decays to 
pK~ir + (2.8%) and semi-leptonically with a 4.5% branching ratio. The heavier ground state 
charm baryons and their excited states (S c and higher) decay via the A c s. The charm-strange 
baryons are assumed to be a negligible contribution to the total cross section. 

Extracting the total charm cross section from data is a non-trivial task. To go from a finite 
number of measured D mesons in a particular decay channel to the total cc cross section one 
must: divide by the branching ratio for that channel; correct for the luminosity, ajj — Njj/Ct; 
extrapolate to full phase space from the finite detector acceptance; divide by two to get the 
pair cross section from the single Ds; and multiply by a correction factor [3] to account for the 
unmeasured charm hadrons. There are assumptions all along the way. The most important is 
the extrapolation to full phase space. Before QCD calculations were available, the data were ex- 
trapolated assuming a power law for the x f distribution, related to the longitudinal momentum 
of the charm hadron by xp = p z /(\ /r S/2) — 2jtit svahy / y/S . The canonical parametrization 
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(1 — xf) c was used where c was either fit to data over a finite xp range or simply assumed. 
These parametrizations led to large overestimates of the total cross section when < c < 2 was 
assumed, especially when data were taken only near xf — 0. Lepton measurements were more 
conservative but were typically at more forward xp. 

Rather than assess the uncertainties in the data, here we address the theoretical uncertainties 
in the calculation of the charm and bottom cross sections. Since the data are generally taken 
in a finite kinematic region, we begin with the calculation of the inclusive distributions to the 
Fixed-Order Next-to-Leading Logarithm (FONLL) level and then discuss the total cross section 
calculations to next-to- leading order (NLO), the most accurate calculation of the total cross 
section over all energies. 

We calculate the transverse momentum, px, distributions of charm and bottom quarks, the 
charm and bottom hadron distributions resulting from fragmentation and, finally, the electrons 
produced in semi-leptonic decays of the hadrons [4]. We then calculate the total charm and 
bottom cross sections, both by the integral over the inclusive px distribution and by integrating 
the total partonic cross section. At each step, we clarify the theoretical framework as well as 
the parameters and phenomenological inputs. Our final prediction is thus not a single number 
but rather an uncertainty band which has a reasonably large probability of containing the 
'true' theoretical prediction. We show that applying this procedure blindly may lead to an 
apparent discrepancy in the two methods, particularly for charm production. We explain why 
this seems to be the case as well as why, when the calculations are done consistently, there is no 
discrepancy. The theoretical uncertainties in both methods of obtaining the total cross section 
are estimated as extensively as possible. We show that, for charm production, the theoretical 
uncertainty on the total cross section is difficult to quantify in a reliable way. 



2 Total heavy flavor cross sections from integrated inclusive distributions 

We first discuss how the total cross section and its accompanying uncertainty is obtained from 
inclusive single particle distributions. We begin with a description of the calculated single 
electron spectrum since heavy flavor hadrons are often observed through their semi-leptonic 
decays, particularly at colliders where direct reconstruction of heavy flavored hadrons at low 
px is difficult. Reconstructed D and B meson decays can only be used to obtain the total 
heavy flavor cross section if they are measured down to px — 0, difficult at colliders. Since the 
RHIC experiments are designed to measure low px hadrons, the full px distribution can be 
accessed. Thus STAR has reconstructed D° — > K^tt? decays to px ~ in addition to their 
single electron measurement [l] . PHEN1X has measured the single electron spectra from heavy 
flavor decays alone [2]. 

The theoretical prediction of the electron spectrum includes three main components: the px 
and rapidity distributions of the heavy quark, Q, in pp collisions at y/S — 200 GeV, calculated 
in perturbativc QCD; fragmentation of the heavy quarks into heavy hadrons, Hq, described by 
phenomenological input extracted from e + e~ data; and the decay of Hq into electrons according 
to spectra available from other measurements. This cross section is schematically written as 

^ = ^^*D{Q^H Q )Vf{H Q ^e), (1) 

where ® denotes a generic convolution. The electron decay spectrum term, /(Hq — > e), also 
implicitly accounts for the proper branching ratio to leptons. 

The distribution EcPctq/cIpq is evaluated at the FONLL level, implemented in Ref. [5]. 
In addition to including the full fixed-order NLO result [617] . the FONLL calculation also 
resums [8] large perturbative terms proportional to a™ log fc (px / m) to all orders with next- to- 
leading logarithmic (NLL) accuracy (i.e. k = n, n — i) where m is the heavy quark mass. 

The perturbative parameters are the heavy quark mass and the value of the strong coupling, 
a s , while the parton densities are a nonperturbative input. We take central values of 1.5 GeV 
for charm and 4.75 GeV for bottom and vary the masses between 1.3 and 1.7 GeV for charm 
and 4.5 and 5 GeV for bottom to estimate the resulting mass uncertainties. 
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Since the FONLL calculation treats the heavy quark as an active light flavor at px >> m, 
the number of light flavors used to calculate a s includes the heavy quark, i.e. n\f + 1 where, for 
charm, n\{ — 3 (u, d and s). The same number of flavors, n\f + 1, is also used in the fixed-order 
scheme where the quark mass is finite. However, in other fixed-order calculations, e.g. to leading 
and next-to- leading order, the number of light flavors is fixed to n\{. The QCD scale at five 
flavors, yl (5) , is set to 0.226 GeV, as in the CTEQ6M parton densities [9]. 

The perturbative calculation also depends on the unphysical factorization (/if) and renor- 
malization (/ij?) scales. The sensitivity of the cross section to their variation can be used to 
estimate the perturbative uncertainty due to the absence of higher orders. We have taken 
Hr,f = Mo = \Jv\- + m2 as a central value in the inclusive distributions and varied the two 
scales independently within a 'fiducial' region defined by /ir^f — £,r,f/J-o with 0.5 < £r,f < 2 
and 0.5 < S,r/S,f < 2. In practice, we use the following seven sets: {(£r,£f)} = {(1,1), (2,2), 
(0.5,0.5), (1,0.5), (2,1), (0.5,1), (1,2)}. The uncertainties stemming from mass and scale varia- 
tions are added in quadrature. The envelope containing the resulting curves defines the uncer- 
tainty. 

The fragmentation functions, D(c — > D) and D(b — > B), where D and B indicate a generic 
admixture of charm and bottom hadrons, are consistently extracted from e + e~ data in the 
context of FONLL [TO] . Using the Peterson et al. fragmentation function [TTJ , with standard 
parameter choices e c ~ 0.06 ± 0.03 and ~ 0.006 ± 0.003, does not provide a valid description 
of fragmentation to FONLL. 

The measured spectra for primary B — > e and D — > e decays are modeled and assumed 
to be equal for all bottom and charm hadrons respectively. The contribution of electrons from 
secondary B decays, B — > D — > e, was obtained by convoluting the D — > e spectrum with a 
parton- model prediction of b — » c decay. The resulting electron spectrum is very soft, giving 
a negligible contribution to the total. The decay spectra are normalized using the branching- 
ratios for bottom and charm hadron mixtures 12J: BK(B — > e) = 10.86 ±0.35%, BR(D — > e) = 
10.3 ± 1.2%, and BR(B -> D -> e) = 9.6 ± 0.6%. 

We first present the transverse momentum distributions for charm and bottom quarks on 
the left-hand side of Fig.[TJ The theoretical uncertainty bands for the two distributions in Fig.[TJ 
is obtained by summing the mass and scale uncertainties in quadrature so that 



where C is the distribution for the central value, /i, max (/i, min) is the maximum (minimum) 
cross section obtained by choosing the central value with the scale factors in our seven fiducial 
sets, and to, max (to, min) is the maximum (minimum) cross section obtained with £r — = 1 
and the lower and upper limits on the quark mass respectively. There is, however, considerable 
arbitrariness in the choice of the method used to assess the theoretical uncertainties. In fact, 
the meaning of the theoretical error due to unknown higher order effects is, to a large extent, 
subjective. The recipe we follow is often used in calculations of cross sections at hadron colliders 
and is similar to the one used to compute heavy flavor cross sections at the Tevatron (see 
Refs. |13|14|15j ). By experience, we assign a probability of 80-90% that the true result lies 
within the band. 

Note that the charm quark uncertainty band is enlarged at low px due to the large value of 
a s at low scales and the increased sensitivity of the cross section to the charm quark mass. In 
Ref. [4], we also noted that, due to the fairly hard fragmentation function, the D meson and 
c quark distributions begin to differ outside the uncertainty bands only for px > 9 GeV while 
the b quark and B meson bands overlap over all px ■ 

The single electron uncertainty bands from D — > e, B — > e and B — > D — > e decays as well as 
the sum are compared to the STAR [TJ and PHENIX data on the right-hand side of Fig. [TJ As 
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Fig. 1. Left-hand side: The FONLL pr bands for c and b quark production compared to the STAR 
D meson data pQ. Right-hand side: The heavy flavor contributions to the single electron spectra: 
D — > eX (dotted), B — ► IX (dashed), B — > D — > IX (dot-dashed) and the sum (solid), is compared to 
the PHENIX [2] and STAR [I] data. 



expected, B —> D — » e is a negligible contribution to the total. While D — > e decays dominate 
at low pt, the B — > e contribution begins to dominate at higher px- The two uncertainty bands 
cross each other in the region 3.5 < px < 12 GeV. The region of crossover is rather broad 
since we consider the c and b quark mass and scale uncertainties to be uncorrelated. If the scale 
uncertainties were assumed to be correlated, the crossover region would be narrower, as shown 
in Ref. |16j . However, for a true measure of the uncertainty, we cannot assume that the scales 
are correlated. The PHENIX measurement is in relatively good agreement with the upper edge 
of the uncertainty band in Fig. [1] while the STAR data tend to lie a factor of 4-5 above the 
central value, falling well above the band. 

If the distributions shown here are integrated over all phase space, the 'perturbative' inputs 
used in the calculation lead to a FONLL total cc cross section in pp collisions of 

^FONLL = 25 g+400 ^ (4) 

at y/S = 200 GeV [4] . The corresponding NLO prediction [4] is 

^NLO = 244 +381 ^ (5) 

The theoretical uncertainty is evaluated as described above. Thus the two calculations are equiv- 
alent at the total cross section level within the large perturbative uncertainties, as expected. 
The total cross section for bottom production is [3] 



(G) 
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Because the FONLL and NLO distributions tend to coincide at small pr and the total cross 
section is dominated by the low px region, the total cross sections and their uncertainties are 
nearly equal in the FONLL and NLO approaches. Earlier papers [17] used m — 1.2 GeV and 
p b. = Pf = Zy/Px + m2 a s reference parameters for charm production. With this choice we 
find <7^ LO = 427 fib, within the calculated theoretical uncertainty band. 

3 Total heavy flavor cross section from total partonic cross sections 

The total partonic cross section has only been completely calculated to NLO [6] . Some NNLO 
calculations arc available near threshold, applicable for y/S < 20 - 25 GeV pTPQj . The NLO 
corrections to the leading order (LO) cross sections are relatively large, K tb = ctnlo/^lo ~ 
2 — 3, depending on p, m and the parton densities ^ZU\. The NNLO corrections are about as 
large at next-to-next-to-leading logarithm [18] but decrease to less than K t h when subleading 
logs are included |19j . Scaling functions [3] proportional to logs of p 2 /m 2 are used to calculate 
the total cross section to NLO. 

The hadronic cross section in pp collisions can be written as 



where x\ and X2 are the fractional momenta carried by the colliding partons and ff are the 
proton parton densities. The partonic cross section is 



where p = 4m 2 / s and ' are the scaling functions to NLO [5]. 

At small p, the 0(a 2 ) and 0{or s ) qq and the 0(a 2 ) gg scaling functions become small while 
the 0{a\) gg and qg scaling functions plateau at finite values. Thus, at collider energies, the 
total cross sections are primarily dependent on the small x parton densities and phase space. 

The total cross section does not depend on any kinematic variables, only on the quark mass, 
m, and the renormalization and factorization scales with central value pr,f = po = m. The 
heavy quark is always considered massive in the calculation of the total cross section and is 
thus not an active flavor in the production calculation. Therefore, the number of light quark 
flavors, uif, does not include the heavy quark while the FONLL calculation uses mt + 1 flavors 
since the heavy quark is an active flavor at high px, as described in Section O 

The theoretical uncertainty on the total cross section is studied within the same fiducial re- 
gion as the pt distributions with the upper and lower limits of the uncertainty band determined 
as in Eqs. j2]) and ([3]). The energy dependence of the charm and bottom total cross sections is 
shown in Figs. [2] and [3] respectively. The left-hand sides of the figures blow up the fixed-target 
and CERN ISR energy regime where the most data are available while the right-hand sides 
show the extrapolation of the cross sections to the collider regime. Only a subset of the most 
recent fixed-target charm data are shown on the right-hand side of Fig. [2] The central value 
of the band is indicated by the solid curve while the upper and lower edges of the band are 
given by the dashed curves. The dotted curves in Fig. [2] are calculated with pp = PR = 2m and 
m = 1.2 GeV, used in Ref. [17] ■ Note that the charm uncertainty band broadens as the energy 
increases. The lower edge of the charm band grows more slowly with y/S above RHIC energies 
while the upper edge is compatible with the reported total cross sections at RHIC [112] . 



With nif light flavors and a fixed scale, the charm and bottom NLO total cross sections at 
y/S = 200 GeV are 





(8) 



NLO 



if 
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(10) 
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Fig. 2. The NLO total cc cross sections as a function of y/S for y/S < 70 GeV (left-hand side) and 
up to 14 TeV (right-hand side) calculated with the CTEQ6M parton densities. The solid curve is the 
central result; the upper and lower dashed curves are the upper and lower edges of the uncertainty 
band. The dotted curves are calculations with m = 1.2 GeV, fiF = = 2m. 




Vs" [GeV] Vs" [GeV] 

Fig. 3. The NLO total bb cross sections as a function of yS for yS < 70 GeV (left-hand side) and 
up to 14 TeV (right-hand side) calculated with the CTEQ6M parton densities. The solid curve is the 
central result; the upper and lower dashed curves are the upper and lower edges of the uncertainty 
band. 

While the central values are only about 25% and 10% higher respectively than the FONLL 
results in Eqs. (j4]) and ([6]), the uncertainty is considerably larger, especially for charm. We now 
discuss the major sources of the theoretical uncertainty and how the apparent discrepancy in 
the magnitude of the charm cross section uncertainty in the RHIC results comes about. 



4 Comparison and discussion 

From the results in the previous two sections, it seems that the total cross section is different 
depending on whether it is calculated from the integral over the inclusive pr distribution or from 
the total partonic cross sections. The difference seems especially large for charm production. 
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Table 1. The values of a a for charm and bottom production at the given values of £r — fin/in. 

£r n u = 3, m = 1.5 GeV n lt = 4, m = 4.75 GeV 
0.5 0.6688 0.2822 

1 0.3527 0.2166 

2 0.2547 0.1804 



This is largely due to the way the strong coupling constant is calculated and the low x, low 
scale behavior of the parton densities. 

In this section, we discuss these two contributions to the theoretical uncertainty and show 
that, if the total cross section is calculated the same way, the two results are, in fact, equivalent, 
as they should be. 



4.1 Strong coupling constant dependence 

The most trivial difference in the two calculations is that the pt distribution is calculated with 
a running scale proportional to ttit while the total cross section is calculated with a fixed scale 
proportional to to. The charm quark uncertainty band is wider at low pr, as shown in Fig. [TJ 
because px < to and the calculation is more sensitive to the lower scale in a s since tot ~ to at 
low pt- While it is more appropriate to use the running scale to calculate inclusive distributions, 
the difference between a fixed and a running scale can be checked by fixing the scale in the 
Pt distributions. The integral of the inclusive distribution increases about 20% for charm and 
about 10% for bottom when a fixed scale is used. This difference is approximately large enough 
to account for the difference in the central values of the total cross section. 

One obviously important contribution to the uncertainty is the difference in the number of 
flavors in the two calculations, especially for charm since the fiducial range, 0.5 < £r < 2, is in 
a region where a s is changing rapidly with hr. Although increasing the number of light flavors 
involves more than just changing a parameter in the calculation of a s , we can get an estimate of 
the importance of the value of a s to the uncertainty in the total cross section by looking at the 
dependence of a s on the renormalization scale. When calculated with the 5 flavor QCD scale for 
CTEQ6M, vlg = 0.226 MeV, and using a scheme where a s is continuous across mass thresholds, 
we have the values shown in Tabled] It is clear, based on these values alone, that the charm 
uncertainty is larger than that for bottom since a s (^u = 0.b)/a s (^R — 2) — 2.63 for charm 
and 1.56 for bottom. The real difference in coupling strength between the two heavy quarks 
is even larger since the leading order cross section is proportional to cq while the next-order 
contribution is proportional to a?. 

Using n\[ + 1 in the FONLL and NLO calculations of the inclusive distributions in Section [2] 
reduces the uncertainty. When the total cross sections in Eqs. (U]) and ([5]) are instead calculated 
with nif , the uncertainty is increased so that the upper and lower limits of the charm uncertainty 
are in agreement with Eq. (J9j) [2.1]. Thus whether charm is treated as a heavy (nif) or an active 
(n\f + 1) flavor in the calculation turns out to be one of the most important influences on the 
limits of the charm uncertainty. 



4.2 Parton density dependence 

Next, we discuss the influence of the parton densities on the theoretical uncertainty. Since to is 
the only perturbative scale, the total cross section calculations in Section [3] are more sensitive 
to the low x and low /i behavior of the parton densities. Probing the full fiducial range of the 
uncertainty band is problematic for charm production since £p = 0.5 is below the minimum 
scale of the CTEQ6M parton densities, /i TE< ^ 6M = 1.3 GeV. Thus, for this scale, backward 
evolution of the parton densities is required. 

The CTEQ6M (NLO, MS scheme) gluon distributions in the fiducial region of the factoriza- 
tion scale, 0.5 < £p < 2, are shown in Fig. 0] The behavior of the gluon distributions for charm 
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Fig. 4. The CTEQ6M parton densities as a function of x for £ = 0.5 (dot-dashed), £ = 1 (dashed) and 
£ = 2 (solid) for m — 1.5 GeV (left-hand side) and 4.75 GeV (right-hand side). The vertical line is the 
value x = 2m/VS in VS = 200 GeV pp collisions at RHIC. 

(left) and bottom (right) are quite different. Since the range 0.5m < fip < 2m for bottom 
quarks lies well above /i TEC ^ 6M , the scale dependence of the gluon density is typical. The gluon 
density increases with decreasing x. The highest low x gluon density is at the largest scale. For 
x values larger than that of central rapidity at RHIC, to the right of the vertical lines in Fig. 2J 
the gluon densities are rather similar although the density is larger at the lower scale. This x 
dependence is quite typical for large, perturbative factorization scales and demonstrates why 
the bottom quark cross section is well behaved as a function of yf~S. 

The gluon distributions with £p = 2 for charm and £p = 0.5 for bottom are similar because 
Hf = 3 GeV and 2.375 GeV respectively. Thus at the highest £p for charm, the low x gluon 
density is well behaved. However, the behavior at lower scales is quite different, especially for 
x < 10~ 2 . When £p = 1, the dashed curve on the left-hand side of Fig. 0] is no longer increasing 
with decreasing x but, instead, is almost flat for x < 10~ 2 with a slight dip in the middle. Lower 
x values are not shown for £p — 0.5 because the backwards evolution gives xg(x,^F = 1) = 0, 
accounting for the high yf~S behavior of the lower bound on the uncertainty band. The low x, 
low /ip behavior of the gluon density depends strongly on how the group performing the global 
analysis chooses to extrapolate to unmeasured regions. All that is required is minimization of 
the global x 2 and momentum conservation. 

4.3 Scale dependence 

Finally, we describe the scale dependence of the charm and bottom cross sections in some detail. 
4.3.1 Bottom 

We first focus on bottom production because the factorization scale is larger than fi Q TE Q 6M j n 
the entire fiducial region. In general, when the factorization scale in Eq. ([7]) is sufficiently high, 
the lowest scales give the highest cross sections because the evolution to higher scales reduces 
the gluon density at x > 10~ 2 . The largest cross sections in the fixed-target regime are thus 
obtained with the combinations (£r,£f) = (0.5,0.5), (0.5,1) and (1,0.5). At relatively low VS, 
the slightly higher value of a s (£,R = 0.5) compensates for the lower gluon density with £f — 1 at 
large x so that the cross section with (£r,£f) = (0-5, 1) is higher than that with (1,0.5). As the 
energy increases and lower x values are probed, the lower small x gluon density with £p = 0.5 
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can no longer overcome the difference between a s (^R = 1) and a s (£,R = 0.5) and, at LHC 
energies, the cross section with (£r,£f) = (1,0.5) is smaller. Indeed, the lower factorization 
scale with (£r,£f) = (0.5,0.5) is not enough to keep the cross section with this parameter 
set larger than that with (1,0.5) or, for that matter, those with either or both £r,£f = 2 at 
sufficiently low x. 

A similar effect occurs for (£r,£f) = (2, 2), (2,1) and (1,2) except that now, at fixed-target 
energies, the set £f) = (2, 2) gives the lowest cross section of the three pairs. In this case, 
(£r,£f) = (1)2) gives the highest cross section of the three sets since the gluon density with 
£f — 1 is higher at large x (lower yf~S) and o, s (£r — 1) is larger, compensating for the slightly 
lower gluon density at large x. At collider energies, £f) = (1, 2) still gives the largest cross 
section of these three parameter sets since the evolution at low x (large y/S) is the dominant 
behavior. At large y/S, the cross section with (£r,£f) = (2,1) drops below those calculated 
with the other two sets. 

These subtle changes in which (£r, £f) set dominates the upper and lower limits of the 
bottom quark total cross section uncertainty band as a function of y/S do not significantly 
broaden the uncertainty band, even at the highest energies because (If > /!q TE( ^ 6M . 

The scale choice in the parton densities affects the dominance of a particular parameter set 
(£kj €f) i n the pt distributions to a lesser extent because at px > m the scales are all large and 
perturbative. Different parameter sets dominate the pt distribution because, at RHIC, higher 
probes the large x range of the gluon distribution while for pt — > 0, x is relatively small. At 
Pt — > 0, the upper and lower edges of the band are thus determined by (£r, £f) = (0.5, 1) and 
(1,0.5) respectively, as is also the case for the total cross sections. However, as px increases, the 
upper and lower edges of the band are defined by (£k,£f) = (0.5,0.5) and (2,2) respectively. 
Increasing pt has the same effect as moving to smaller y/S : both probe larger x where the gluon 
distribution with £f = 0.5 is higher than that with £f = 2, as is obvious from the right-hand 
side of Fig. [4j 

Even though the scale dependence of bottom production is not negligible, as we have seen, 
it is not strong in the defined fiducial range. The difference in the bb cross sections in Eqs. ([6]) 
and (|10p can be almost entirely attributed to the change from the running scale in Eq. |T]) and 
the fixed scale in Eq. |(7J| . Thus the bottom production cross section is rather well under control. 



4.3.2 Charm 



However, the scale dependence of the total charm cross section on y/S is another story due to 
the behavior of the CTEQ6M gluon distribution at charm quark scales. Since m ~ /j,g TE( ^ 6M , 
using the full fiducial region to estimate the theoretical uncertainty on the total charm cross 
section problematic. The smaller charm mass exaggerates the factorization scale dependence of 
the total cross section described above for bottom production. 

Thus the charm quark uncertainty band on the total cross section, Fig. [2 spans an order of 
magnitude at fixed-target energies, increasing to the value given in Eq. ([9]) for nif at y/S = 200 
GeV. The low scale behavior for (£r,£f) = (0-5,1) and (1,0.5) defines the upper and lower 
edges respectively of the uncertainty band at collider energies. Indeed, for the total cross section 
calculated with n\{ light quark flavors, the STAR point [1] is compatible with the upper limit 
of the band although the inclusive px data lies above the FONLL calculation with n\f + 1 light 
flavors [22] • However, we stress that this apparent agreement of the STAR result with the total 
cross section does not mean that the discrepancy between the high px STAR results and the 
FONLL prediction can be ignored. At high px, the FONLL calculation is more reliable since 
here charm is correctly treated as an active flavor, with nu + 1, and light quark effects are 
resummed, improving the prediction at finite px- 

The charm band grows broader with increasing y/~S, corresponding to decreasing x. At 10 
TeV, the width of the uncertainty band has increased to almost two orders of magnitude. Thus, 
without a better handle on the gluon density at low x and low scales, one may question whether 
such a large uncertainty is meaningful. It may also be questionable whether the lowest scales, 
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£-R> £f = 0-5 should be included in the calculation of the charm uncertainty, especially when 
Hf < /x T 6M for three light flavors. 

The full fixed-target data set also exhibits a large uncertainty due to the method of extrap- 
olation used, the assumed branching ratios and the A dependence, as shown on the left-hand 
side of Fig. [21 However, if only the most recent data are used, the uncertainty in the data seems 
to be reduced. As an alternative, one may try to 'fit' the mass and scale parameters to these 
data [17] for /i > m. The dotted curves in Fig. [2] show the energy dependence of one such 
attempt with m = 1.2 GeV, {£r,£,f) = (2,2). The calculated cross section lies just above the 
central value of the band and, although the quark mass is smaller than the assumed central 
mass value, the larger value of £p guarantees a more regular y/S dependence than that obtained 
with smaller values, as shown in Fig.UJ 

5 Conclusions 

We have shown that when the total cross section is calculated with the same parameter sets 
and the same number of light quark flavors, a consistent result is obtained by both integrating 
over an inclusive distribution and starting from the total partonic cross section, as should be 
expected. However, the charm results are extremely sensitive to the number of flavors, the scale 
choice and the parton densities. One of the biggest sources of uncertainty in the total charm 
cross section at collider energies is the behavior of the gluon density at low x and low scale, as 
yet not well determined. Until it is further under control, better limits on the charm quark total 
cross section will be difficult to set. A complete NNLO evaluation of the total cross section may 
reduce the scale dependence but will still be subject to the same types of uncertainties. 

It is thus not clear which estimate of the total charm cross section uncertainty, Eqs. (JU) 
and (jSJ) or Eq. ([9]), is more reliable. If the low px region is ignored and the heavy quark may 
be considered an active flavor then the appropriate number of flavors is Tin + 1 rather than n\i 
and the smaller error band used to compare the RHIC pr distributions 22J is more reasonable. 
However, when heavy flavor production is measured over the full pr range, down to pt ~ 0, 
then three light flavors should likely be used for charm, resulting in the larger uncertainty. 
Unfortunately, in this case, the uncertainty is driven by scales lower than the initial scale of the 
parton density, further complicating the interpretation of the limits on the uncertainty band. 
Thus, rather than arbitrarily choosing one result over another, we prefer to stress that there is 
little predictive power in the charm production uncertainty. 
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